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Introduction 


The  goal  of  the  effort  reported  here  is  to  provide  a  convenient,  efficient, 
and  accurate  method  to  estimate  the  rate  of  single  event  upsets  for  devices  in 
the  space  environment  coring  the  various  phases  of  the  CRRES  mission.  This 
calculation  integrates  results  from  the  CRRES  cosmic  ray  model,  results  from 
ground  tests  of  device  upset  susceptibility,  approximate  device  geometry 
moctels,  and  device  shielding  estimates  to  obtain  the  upset  rate.  The  method 
provides  upper  and  lower  limK  estimates  of  the  upset  rate  from  different 
integrations  of  the  upset  cross-section  test  results.  These  limits  provide  some 
incfication  of  the  precision  of  the  result  and  its  sensitivity  to  the  quality  of  the 
device  test  dsrta. 

This  report  provides  an  overall  description  of  the  approach,  a  detailed 
description  of  the  specific  method  adopted,  and  full  documentation  of  the 
software  implementation  of  that  method. 


Methods  and  Procedures 

The  method  chosen  to  calculate  the  SEU  rate  is  similar  to  that  described 
by  Pickel  and  Blandford  (1980).  It  is  a  direct  application  of  a  few  simple 
geometrical  concepts. 

The  upset  rate  is  determined  as  the  rate  at  which  energetic  heavy  ions  in 
the  space  radiation  environment  can  deposit  more  than  a  certain  minimum 
amount  of  energy  into  a  specific  volume.  The  minimum  energy  deposit  required 
to  cause  an  upset  is  called  the  upset  threshold.  The  specific  volume  is  the 
sensitive  volume  of  a  cell  in  the  ^vice.  The  particle  rate  is  determined  from  the 
procfcjct  of  the  local  particle  flux  (per  unit  time  per  unit  solid  angle  and  per  unit 
area)  and  the  geometrical  factor  presented  by  the  sensitive  volume  (solid  angle, 
area  product). 

The  requirement  to  deposit  more  than  a  certain  minimum  amount  of 
energy  into  a  specific  volume  restricts  the  c^culation.  It  leads  to  a  specific 
method  to  organize  the  particle  flux  distribution  and  a  specific  method  to 
describe  the  geometry  of  the  sensitive  region.  This  modified  description  of  the 
particle  flux  distribution  is  called  the  linear  energy  transfer  (LET)  spectrum.  The 
modified  description  of  the  geometry  of  the  sensitive  region  is  called  the  path 
length  distribution. 

LET  ^MCtrum  The  linear  energy  transfer  (LET)  of  an  energetic  heavy  ion,  the 
energy  loss  per  unit  distance  travelled  in  a  material,  generally  increases  with 
the  atomic  number  of  the  ion  (Z)  and  generally  decreases  wKh  increasing 
velocity  (except  at  low  and  high  velocities  where  atomic  effects  and  relativistic 
effects,  respectively,  become  important).  'Hius  a  fast  iron  ion  and  a  slower 
oxys^n  ion  can  have  the  same  LET  and,  by  assumption,  the  same  ability  to 
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produo*  •  event  In  a  specific  device.  The  LET  spectoim  is  a 
convenieftt  way  fi)  inteqwle  and  track  of  the  contributiorrs  of  various  ions 
ORJeied  ttcooRinQ  to  LET.  The  integral  LET  spectrum  used  here  is  reported  as  a 
product  of  the  CRRES  Cosmic  Ray  Model  Group,  it  represents  the  total  flux  of 
pardciee  of  types  «ftth  an  LET  larger  than  wme  threshold  (L).  This  is  written 
a8F(D. 

PbNi  length  dietrNMitkMi  The  energy  deposited  by  a  partide  with  LET  Lo 
along  path  I  Is  ^  »  Lo*i.  Thus  the  other  information  needed  to  establish  the 
energy  deported  an  ion  penetrating  a  region  is  the  length  of  the  chord 
through  the  region  0)  r^Mesenting  the  track  of  the  ion. 

Viewed  horn  an  Mbhrwy  rSrection,  any  solid  presents  a  distribution  of 
chord  lengtfw  or  projected  thicknesses,  over  its  projected  area.  This  is  the 
unidirectional  pa^  length  distribution  through  the  solid.  It  has  a  simple 
interpretttion  but  can  be  rather  complicated  to  express  mathematically. 

For  any  solid  at  an  arbfirary  orientation,  there  is  an  area,  solid  angle 
product  (georncrtricel  factvV)  presented  by  the  solid  for  a  given  path  length.  In  the 
present  work,  since  the  region  of  interest  has  an  arbitrary  orientation  in  space, 
the  path  length  distrttxjtion  is  computed  as  the  average  over  all  orientations. 

The  path  ler^  cRstrttHJtion  is  wrMen  P(i).  It  is  zero  for  path  lengths  longer  than 
the  maximum  chord  thrcHjgh  the  volume  (Imax)*  For  the  work  reported  here,  the 
dshttMJtion  is  normaHzed  so  that  the  integral  over  all  path  lengths  is  equal  to  the 
total  geometricai  factor  of  the  solid.  For  many  types  of  solid  prisms  this  is  x  times 
the  surface  area  of  each  face  of  the  soNd.  For  a  right  rectangular  solid  prism  with 
dimensions  x.y.2  the  path  length  cNstribution  is  written  as  P(l:x,y,z)  and  the 
normaiizsMon  condition  is: 


P(l,'X,y,z)  cH  »  2jc(xy  +  xz  ♦  yz) 


(1) 


The  maximum  path  length  is  given  by: 

lmax»Vx2  +  y2  +  22  (2) 

The  form  of  the  function  P()  is  provided  by  Rckel  and  Blandford  (1960). 
afihough  tl^  choose  to  use  a  different  normailzation  than  presented  here. 


Device  geometry  Single-event  upset  susceptibility  tests  of  microelectronic 
devices  are  often  characterized  by  two  parameters:  a  cross-sectionai  area  of  the 
device  which  is  sensitive  to  upset,  and  a  threshold  value  of  the  LET  necessary 
to  induce  an  upset  (ref.)  in  the  simplest  Idealization  of  the  test  results,  the  upset 
cross-section  is  zero  below  some  threshold  LET  value,  Lc.  and  a  constant  equal 
to  some  area.  Au.  above  this  threshold. 

The  u|^  cross-section  result^  yield  to  a  simple  geometric  interpretation 
which  has  been  validated  by  continuing  device  tests  and  analyses.  The 
sensitive  region  of  any  celt,  or  bit  of  RAM,  is  often  modelled  as  a  rectangular 
soSd  section  of  ttie  stflcon  chip.  This  section  has  a  thickness  (z)  and  some 


2 


surface  area  (x*y).  The  geometric  interpretation  of  the  total  sensitive  area 
measured  in  the  SEU  tests  is  clear:  For  a  device  with  N  cells  (e.g.  N  bits  in  a 
fWMi  the  total  sen^tive  area  of  the  device  is  Au  =  N*x*y.  Thus,  a  measurement 
of  the  upset  cross  section  together  with  a  knowledge  of  the  number  of  bits  (N) 
yields  an  estimate  of  the  sensitive  cell  area  (x*y)-  The  individual  dimensions  (x 
arid  y)  can  often  be  obtained  with  sufficient  accur^  just  as  the  square  root  of 
that  area.  More  precise  estimates  are  rarety  required  or  justified.  They  must 
come  from  details  of  the  manufacture  of  the  device  or  from  physical  analysis  of 
tH0  struclLH'0 

The  upset  threshold  value  of  the  LET  (Lc)  multiplied  by  the  thickness  of 
the  region  (z)  is  interpreted  as  a  minimum  energy  or  charge  deposit  necessary 
to  change  the  state  of  the  device,  Ec  -  z*Lc.  (For  silicon,  this  energy  deposit  or 
charge  deposit  are  proportional  to  each  other.  ITie  constant  of  pro^rtionaiity  is 
3.6  eV/electron.)  The  thickness  (z)  is  not  generally  known  with  great  precision.  It 
depends  on  parameters  associated  with  the  device  technology  and  wafer 
processing.  Fortunately,  since  both  the  test  results  and  the  real  spectrum  are 
expressed  in  terms  of  the  l£T.  an  impredse  knowledge  of  the  depth  (z)  has  little 
effect  on  the  precision  of  the  result.  This  has  been  demonstrated  by  detailed 
calculation  (Shoga.  et  al.,  1988). 

Total  upset  rate  The  total  upset  rate  (U)  for  the  simple  model  outlined  above 
is  given  by  the  following  expression: 


U-  F(Ec/l)P(l:x.y.z)dl 

Jo  (3) 

The  integral  is  performed  over  the  entire  path  length  distribution.  However,  for 
any  path  length,  I,  only  particles  with  LET  above  the  threshold  L  =  Ec  / 1 
contribute  to  the  upset  rate.  Because  the  LET  spc^rum  and  the  path  length 
distribution  are  not  simple  analytic  functions,  the  integral  is  generally  computed 
numerically. 


Integrating  Detailed  Upset  Cross-Sections 

The  relatively  simple  integral  described  in  the  previous  section  is 
apprc^ate  in  the  i^aiized  situation  where  the  measured  upset  cross  section  is 
a  step  function,  equal  to  zero  below  some  threshold  value  of  the  LET  and  equal 
to  a  constant  value  above  it.  The  situation  becomes  slightly  more  compHcated  in 
the  real  world,  in  this  section  other  shapes  for  the  upset  cross-section  are 
addressed. 

The  upset  cross-section  measurements  are  rarely  step-functions.  More 
often,  there  is  a  transition  region  of  LET  where  the  cross-section  rises  from  zero 
to  some  asymptotic  value  (Figure  1).  This  behavior  has  been  discussed 
extensively  in  the  literature  and  various  explanations  offered.  For  our  purposes, 
however,  we  choose  to  actopt  another  simple  concept  to  understand  this  effect. 
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Tito  con(to>l  diso  motivates  and  justifies  our  method  to  incorporate  the  effect 
into  the  calcutedion. 


Figure  1 :  Measured  single  event  upset  cross  sections  for  the  93L422  RAM 


The  concept  to  apply  to  the  variation  in  upset  cross-section  with  LET  is  to 
treat  the  ratio  of  the  upset  cross-section  to  its  asymptotic  value  as  an  upset 
“effectiveness”  or  “efficiency”.  To  accomplish  this  the  upset  cross  section  is 
divided  by  the  area  factor  (x*y).  At  low  values  of  the  LET  the  efficiency  for  upset 
in  the  sensitive  region  is  low.  As  the  LET  increases,  the  efficiency  also 
increases.  As  the  LET  is  increased  still  more  an  asymptotic  cross  section  is 
established  which  defines  the  area  sensitive  to  the  SEU  process.  The  upset 
cross-section  divided  by  its  asymptotic  value  is  an  efficiency  in  the  range  from  0 
to  1.  The  area  factor  (x*y)  is  retained  through  the  path  length  distribution 
normaHzation  as  in  the  discussion  above. 

In  applying  this  concept  to  the  integral  necessary  to  calculate  the  SEU 
rate,  the  cross-section,  renormalized  as  an  efficiency,  can  be  treated  as  a  sum 
of  step  Unctions  in  the  LET  (h(L))  with  amplitude  Sj  at  location  Lj: 

e(L)  =  ^-5;iSih(Li)-Si.,h(L|)l 

^  i  (4) 

The  subbacdon  is  necessary  to  remove  the  contribution  of  one  step  function 
from  the  region  to  the  right  of  the  next.  Alt  of  the  values  Sj  are  in  the  range  from  0 
to  1  arKf  rmmoionicalty  increasing  with  increasing  L|.  The  sequence  of  values  si 
and  Li  can  be  ch(»en  in  many  ways  to  approximate  the  shape  of  the  upset  cross 
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section.  Lower  and  ujii^r  limit  estimates  can  be  chosen  and  refined  at  will.  An 
example  is  illustrated  in  Figure  2. 


Lower  Limit 


Upper  Limit 


i-i-i  Lj  let 

Figure  2:  Step  function  approximations  to  the  upset  cross  section 

With  this  formulation  of  the  measured  upset  cross  section,  the  complete 
forrmila  for  the  upset  rate  becomes: 

U- Y  (Sj-Sj-i)!  F(zLj/l)P(l:x,y,z)dl 

i  Jo  (5) 

The  tector  z*Li  is  the  critical  energy  deposit  threshold  for  each  step  function 
piece.  This  completes  the  presentation  of  the  method. 
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Application  of  the  method 


Software  to  implement  the  method  described  in  the  previous  section  has 
been  developed.  A  detailed  description  of  the  subroutines  and  listings  of  the 
source  code  are  presented  In  the  appendices.  Some  practical  considerations 
for  the  use  of  this  software  are  provided  in  this  section. 

Following  the  discussion  above,  there  are  four  major  components  in  the 
calculation.  E^h  of  these  components  must  be  provided.  Most  of  them  are  input 
parameters  or  input  data  files.  The  major  components  are: 

1 .  An  estimate  of  the  dimensions  of  the  device  cell  which  is  sensitive  to 
upset,  modelled  as  a  rectangular  solid  prism,  and  the  number  of  sensitive 
regions  In  the  device.  The  dimensions  are  labelled  x,  y.  and  z.  By  convention, 
the  z  dimension  is  normal  to  the  face  of  the  chip  die.  The  area  A  »  x*y  is  the 
maximum  total  sensitive  area  of  a  unit  ceil  (the  sensitive  area  corresponding  to 
a  bit  in  a  RAM.  for  example).  N  is  the  number  of  unit  cells  in  the  device.  (The 
number  of  bits,  for  example.)  The  quantity  N*x*y  should  be  equal  to  the 
asymptotic  value  of  the  upset  cross  section  (part  3,  below). 

2.  The  distribution  of  path  lengths  through  the  sensitive  region,  P(l).  The 
model  provided  in  the  code  is  for  that  of  a  solid  rectangular  prism. 

3.  The  upset  cross-section  of  the  device  as  a  function  of  LET,  a(L).  The 
cross-section  is  obtained  from  ground  tests.  Physically,  to  be  consistent  with  the 
assumptions  of  the  model,  the  upset  cross  section  must  be  monotonicaily 
increasing  with  Increasing  LET.  Often  the  test  data  are  not  monotonic  (c.f.  Figure 
2).  This  is  attributed  to  imperfections  in  the  test  results  and  the  scatter  around  a 
smooth  monotonic  curve  can  be  used  as  an  estimate  of  the  standard  deviation 
of  the  test  data  set.  The  data  set  which  is  to  be  used  in  this  software  must  be 
smoothed  (by  hand,  if  necessary)  so  that  it  is  monotonic. 

4.  The  flux  of  particies  in  the  radiation  environment  as  a  function  of  LET, 
F(L).  This  distribution  is  provided  as  one  of  the  products  of  the  CRRES  Cosmic 
Ray  Model,  calculated  from  models  and/or  measurements  of  the  abundances 
and  energy  spectra  of  the  energetic  heavy  ions  in  the  cosmic  rays.  It  is 
calculated  for  a  variety  of  assumptions  for  the  thickness  of  shielding  between 
the  environment  and  the  sensitive  region  of  the  device.  The  LET  spectrum  for 
the  proper  epoch  and  for  the  proper  thickness  of  passive  shielding  must  be 
selected  to  retain  accurate  upset  rate  estimates. 

For  each  run,  the  upset  rate  is  cateulated  three  times,  once  each  for  a 
lower  Nmit.  an  upper  limit,  and  a  central  estimate.  All  of  these  estimates  are 
controlled  through  the  in^kit  cross  section  data  file.  Both  upper  and  lower  limits 
are  defined  with  respect  to  the  LET  and  cross  section  values  provided  in  the  file. 
In  this  calcutation  of  the  SEU  rate,  the  step  functions  used  for  the  lower  and 
upper  Hmit  estimates  are  placed  at  the  values  of  the  LET  which  are  provided  in 
the  input  file.  TTie  lower  limit  model  uses  the  cross  section  listed  in  the  file  at 
each  LET.  This  j^lds  a  lower  limit  estimate  as  indicated  by  inspection  of  Figure 
2.  The  upfMr  limit  model  uses  the  cross  section  provided  in  the  input  file  at  the 
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next  higher  value  of  the  LET.  This  yields  an  upper  limit  estimate,  again,  as  is 
clear  from  insp^ion  of  Figure  2.  The  central  estimate  is  obtained  by  a  20-point 
linear  interpolation  of  the  listed  cross  section  vs  LET  values.  The  effects  of 
various  uncertainties  in  the  measured  cross  sections  can  be  investigated  by 
modifying  the  cross  section  data  file.  Different  runs  can  be  made  with  fewer 
entries  in  the  cross  section  table  or  with  the  additional  cross  section  values 
obtained  by  any  of  a  variety  of  user-generated  interpolation  schemes. 

Additional  details  of  specific  data  formats  and  the  mechanics  of  using  the 
code  are  given  in  the  appendices. 
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Appendix  A 


CRRES  Cosmic  Ray  Model 
Single  Event  Upset  Rate  Calculator 

SEU_CALC  and  related  support  subroutines 


Source  code  organization  and  layout 

All  of  the  source  code  Is  written  in  ANSI  C.  The  source  code  is  contained 
in  three  files  named  SEU_CALC.C.  PATHDIST.C.  and  INP_LIST.C.  All  of  the 
source  code  is  in  the  public  domain.  The  subroutines  contained  in  each  file  are 
listed  below; 

SEU_CALC.C 

main() 

integrandO 

myfluxO 

IgaussO 

doJntO 

PATHDIST.C 

pIdsetO 

pIdtO 

PW() 

INP.LIST.C 

inputJistO 
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Subroutine  descriptions 


A  brief  description  of  the  purpose  served  by  each  subroutine  is  provided 
here  togMher  with  the  relationships  between  them  and  (filing  sequences. 

Subroutine  name:  main() 

Calling  sec^ence:  None,  entry  point  for  SEU.CALC 
^turned  value:  None 

Description:  main()  is  the  entry  point  for  the  SEU_CALC  pad^age. 

It  opens  the  output  file,  handles  the  user  interactions,  writes  the 
documerrtation  of  the  run.  handles  units  conversions,  and  caiis  the 
other  routines  in  order  as  required. 

Subroutine  name:  integrand() 

Calling  seouence:  double  integrand(double  I).  where  I  is  a  path  length 
in  rmorons  through  the  sensitive  volume  of  the  device 
Returned  value:  a  double-precision  value  of  the  integrand  for  path 
length  I  given  the  path  length  dstrtbution.  the  LET  ^)ectrum,  and  a 
critical  charge  threshold  ^ 

Description:  Given  a  path  length  distribution  p(l).  a  critical  charge 

threshold  Q.  and  an  LET  spectrum  F(L).  integrand(i)  returns  the 
number  of  particies  per  second  capable  of  depositing  more  than  Q 
in  path  length  I  as  pOrF(LsQ/l) 

Subroutine  name:  myftuxO 

CalKng  seci^ence:  double  myflux(double  L).  where  L  is  an  LET 
threshold  in  units  of  picocouiombs  per  micron 
Returned  value:  a  double-precision  value  of  the  integral  flux  with  LET 
greater  than  L  F(L) 

Description:  Given  the  tabulation  of  the  LET  spectrum,  myflux(L) 

rehjms  the  interpolated  value  of  the  spectrum  at  LETsL.  If  L<0  or  L 
Is  greater  than  the  largest  value  in  the  tabulation,  the  flux  returned 
issero. 

Subroutine  name:  tgauss() 

CiMng  sequence:  double  lgauss(doubte  a,  double  b) 

Returned  value:  a  double-precision  value  of  the  integral  of  integrand() 
(above)  from  a  to  b 

Description:  Igauss  does  a  ten-point  Gauss-Legendre  integration 

based  on  Press,  stat.  (1988). 

&rt>routine  name:  doJnt() 

CaMng  sequence:  int  doJnt(doubie  Acc,  double  min.  double  max. 
double  *res,  double  *old) 

Reaimed  ^«due:  an  integer  number  of  integrsrtion  siHxiivisions  used 
as  weH  as  the  current  best  (*res)  and  prevkMis  (*old)  integrals 
Desertion:  doJnt()  integrates  the  value  of  integrand()  (s^ve) 

hwn  min  to  m»(.  The  number  of  subdivisions  of  the  int^ration  is 
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increased  until  two  successive  results  (a1  and  a2)  agree  to  within 
less  than  Acc*a2.  The  integral  is  first  performed  with  a  20-point 
trapezddal  mie.  Then  it  is  repeated  using  Qauss-Legendre 
intregrals  over  smaller  and  smaller  divisions  of  the  range  of 
integration  until  convergence  to  the  required  accuracy  is  achieved. 

Si^routine  name:  pldset() 

Calttng  sequence:  double  pldset(double  x.  double  y,  double  z) 

Returned  vahje:  a  double  value  number  equal  to  the  maximum  path 
length  through  the  sensitive  region 

Description:  pldset()  tabulates  the  path  length  distribution  for  a 

right  rectangular  solid  prism  as  calculated  by  pld().  This  is  done  to 
speed  up  the  calculation.  The  path  length  distribution  is  tabulated 
at  1000  points  between  the  minimum  and  maximum.  When  values 
are  required  later  they  can  be  obtcuned  by  simple  interpolation. 

Subroutine  name:  pldt() 

Calling  sequence:  double  pldt(double  s) 

Returned  value:  the  double  value  of  the  path  length  distribution  at 
path  length  s 

Description:  pldt()  interpolates  the  tabulation  made  by  pIdsetO 

Subroutine  name:  pid() 

Calling  sequence:  double  pldt(double  s.  double  x.  double  y.  double  z) 
Returned  value:  the  double  value  of  the  path  length  distribution  at 
path  length  s  for  a  right  rectangular  solid  of  dimension  x.y.z 
Description:  pid()  implements  the  formula  in  the  appendix  of 

Pickei  and  Blandford  (1980).  pld()  uses  a  different  normalization 
than  do  Pickei  and  Blandford,  however.  In  pld()  the  path  length 
distribution  is  normalized  such  that  the  result  returned  is  the 
geometrical  factor  of  the  right  rectangular  solid  prism  for  chords  of 
length  s.  The  unKs  of  this  result  are  area  *  steracfians  with  the  area 
as  the  square  of  the  units  in  which  x,y,z,  and  s  are  presented.  The 
integral  of  pld()  over  all  values  of  s  is  equal  to  the  total  geometry 
factor  of  the  sofid,  which  is  equal  to  pi  times  the  sum  of  the  areas  of 
the  six  faces. 

Subroutine  name:  inputJistO 

Calling  sequence:  int  input Jist(char  ‘file,  char  ‘key,  int  max,  double  ‘cl , 
double  ‘c2,  char  ‘tnie) 

Returned  \«lue:  the  number  of  items  retrieved  from  the  file 
Description:  input Jist()  reads  the  key-based  input  data  files 

needed  by  the  program.  The  file  name  is  contained  in  the  ‘file 
string  and  the  selection  key  in  the  ‘k^  string.  The  ‘title  string  is 
used  to  return  the  comment  following  the  key  (if  maxn  is  non-zero). 
If  maxn  is  zero  then  the  value  returned  in  ‘cl  is  the  numerical 
vakie  immediately  foliowing  the  key.  Otherwise  max  defines  the 
maximum  number  of  pairs  of  entries  to  be  returned  from  the  table, 
cl  and  c2  should  point  to  double  arrays  that  are  at  least  max  long. 
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Tht  values  read  from  the  table  cohimns  in  the  input  file  are 
rMurrted  in  cl .  for  the  first  column,  and  in  c2  for  the  second 
column. 


input  and  output  file  descriptions 

Tmo  fitos  of  input  ctata  are  required  to  perform  the  SEU  rate  caicuiation: 
in  LFT  wedrum  file  (LET_SPECTRUM).  and  a  device  cross-section  file 
(DEVICE).  Both  files  are  text  files  srfth  the  necessary  information  presented 
toNovrlng  key  ttrings  or  in  a  tabular,  columnar  format.  The  contents  and  the 
format  rsquhements  and  restrictions  for  these  files  are  described  here.  In 
addWon  to  writino  to  the  user  terminal  or  console  screen,  the  program  also 
srrtles  its  oufiM  to  a  file  to  provide  a  permanent  record  of  the  mn. 


are  ABCii  files  corrtaining  descriptive  text,  a  key  stHng  of 
chaieoiem.  and  numbers.  Lines  are  delimited  by  either  a  new-Kne  character  or 
by  an  end  of  file.  The  values  may  be  in  any  floating-point  format  recognized  by 
the  ANSI  C  atofO  functton.  The  values  found  in  the  input  files  as  interred  by 
the  pro^am  are  written  to  the  output  file.  When  first  using  a  new  inpd  DEVICE 
or  L^^SPECTRUM  file,  the  output  should  be  carefully  checked  to  ensure  that 
the  interpretation  of  the  data  is  correct  and  complete. 


This  file  must  be  named  or  (in  VMS,  for  exarr^)  defined  or  assigned  to 
the  name  LET_5PEC1THM.  The  key  string  is  the  three-character  strfog  "LET. 
Any  number  of  lines  of  text  may  begin  the  file  as  tong  as  the  sequence  "LET  is 
not  present.  The  first  line  conttiining  "LET  is  assumed  to  be  a  comment  fine  and 
fo  copied  to  me  oumut  as  partial  documentation  of  the  run.  Foiiowing  tNs 
comment  firm  must  be  a  series  of  lines  conta'ning  two  values  per  line  as  in  a 
two-cdumn  tabular  form.  This  series  of  values  is  the  ttfoulation  of  the  LET 
qpsctnim.  The  first  value  on  each  line  must  be  the  LET  in  units  of  MeV/micron. 
1m  second  value  must  be  the  integral  flux  corresponding  to  that  value  of  the 
LET  in  units  of  particles  per  square  meter  per  steradian  per  second.  A  maximum 
of  entries  in  the  LET  spectrum  table  is  supported  at  present.  This  number 
can  be  changed  by  changing  the  value  of  NUMSPC  and  re-compiling.  A 
sample  LET  spectrum  file  is  provided  as  Appendix  B. 


must  be  named  or  (in  VMS,  for  example)  defined  or  assigned  to 
the  name  DEVICE.  There  are  five  sets  of  key  strtngs  which  must  appear  in  the 
fits.  The  first  four  key  strings  are  the  two-character  sequences:  "X«",  *Y«”.  "Z«”, 
and  *N«".  The  iMt  k^  is  the  three-character  string  "LET.  The  "X«",  "Y«".  “ZmT, 
did  "NiB*  key  strings  must  each  be  followed  immediately  (except  for  optional 
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blank  spaces  and  tabs)  by  a  numeric  value.  This  value  is  interpreted  and 
assigned  to  the  respective  variable. 

The  first  three  keys  describe  the  physical  characteristics  of  the  device 
and  Its  sensitive  vokime.  X,  Y.  and  Z  are  ttw  dimen^ons  of  the  sen^tive  region 
ceil  in  microns.  The  Z  dimension  is  assumed  to  correspond  to  depth  in  the 
siHcon.  The  (^entity  X*Y  is  the  sensitive  surface  area  of  the  cell. 

The  number  of  sensitive  volumes  in  the  device  cross-section  table  is 
defined  by  the  value  of  N.  If  the  cross-section  is  tabulated  per  bit.  N  should  be 
set  to  1. 

Any  number  of  lines  of  text  may  begin  the  file  (these  should  include  the 
"Xa”,  ‘Ya”,  "Za”,  and  "Na*  key  strings)  as  long  as  the  sequence  “LET'  is  not 
present.  The  first  line  containing  “LET  is  assumed  to  be  a  comment  line  and  is 
copied  to  the  output  as  partial  documentation  of  the  run.  Following  this  comment 
Kne  must  be  a  series  of  lines  containing  two  values  per  line  as  in  a  two-column 
tabular  form.  This  is  the  tabulation  of  the  device  upset  cross-section.  The  first 
value  on  each  line  must  be  the  LET  in  units  of  MeV  per  milligram  per  square 
centimeter.  The  second  vakie  must  be  the  device  upset  cross-section 
corresponding  to  that  value  of  the  LET  in  units  of  square  centimeters  per  device 
(or  square  centimeters  per  bit  if  the  number  of  bits  is  set  to  1 .  above).  A 
maxiriuim  of  50  entries  in  the  device  upset  cross-section  table  is  supported  at 
present.  This  number  can  be  changed  by  changing  the  value  of  NUMXST  and 
re-oompilir^).  A  sample  device  upset  cross-section  file  is  provided  as  Appendix 
C. 


Output  file  description: 

One  output  file  is  generated  whenever  SEU_CALC  is  run.  It  is  named 
SEU_CALC.  The  output  file  provides  a  version  number  of  the  program,  the 
user^lected  convergence  parameter,  documentation  of  all  of  the  data  used  in 
the  calculation  as  read  from  the  LET_SPECTRUM  and  DEVICE  input  files,  and 
the  results  of  the  integration.  A  sample  output  file  is  provided  as  Appendix  D.  for 
a  run  using  the  input  files  of  Appendices  B  and  C. 
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Appendix  B:  Sample  LET^SPECTRUM  input  file 


LET  (MeV/ micron) 

{  l.OOOOe-02 
{  1.33886-02 
{  1.79256-02 
{  2.39986-02 
{  3.21296-02 
{  4.30156-02 
{  5.75906-02 
{  7.71036-02 
{  1.03236-01 
{  1.38206-01 
{  1.85036-01 
{  2.47736-01 
{  3.31666-01 
{  4.44046-01 
{  5.94496-01 
{  7.95926-01 
(  1.06566-i-OO 
{  1.42676+00 
{  1.91006+00 
{  2.55726+00 
{  3.42376+00 
{  4.58376+00 
{  6.13686+00 
{  8.21616+00 
{  1.10006+01 


Flux  (/m2  s  sr) 

2.19726+01 
1.81846+01 
1.44266+01 
1.11216+01 
8.08896+00 
5.95116+00 
4.38736+00 
3.16576+00 
2.26086+00 
1.64686+00 
1.17566+00 
7.73796-01 
4.0817e-01 
1.99986-01 
1.00806-01 
5.15046-02 
2.66136-02 
1.36886-02 
6.84066-03 
3.4166e-03 
1.53846-03 
6.72406-04 
1.39426-04 
5.6336e-08 
0.0  } 


"LEASAT  1"  sp6ctrum 

}, 

), 

), 

), 

), 

), 

), 

}, 

), 

}, 

}, 

), 

), 

), 

), 

}, 

), 

), 

}, 

}/ 

}, 

), 

), 
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Appendix  C:  Sample  DEVICE  input  file 


This  is  for  the  93L422  RAM 
X-39  Y-39  Z-1 

N-  1024 


/*  Upset  cross-section  data  table. 

The  example  given  here  is  the  measured 
upset  cross  section  for  an  AMD  93L422  RAM 
tested  22  May  1987  at  LBL  88”cyclotron. 
LET  (MeV/mg/cm2)  Cross-Section  (cm2/device)  */ 


{ 

0.1  , 

0.0 

{ 

0.5  , 

l.OOe-4 

I, 

{ 

1.0  , 

3.23e-3 

}/ 

{ 

1.4  , 

4.34e-3 

{ 

3.0  , 

5.77e-3 

{ 

4.2  , 

7,20e-3 

{ 

5.5  , 

8.73e-3 

1/ 

i 

8.2  , 

l.OOe-2 

{ 

9.6  , 

1.12e-2 

{ 

15.  , 

1.18e-2 

}, 

{ 

21.  , 

1.24e-2 

}, 

{ 

30.  , 

1.30e-2 

}, 

{ 

35.  , 

1.39e-2 

}/ 

{ 

42.  , 

1.48e-2 

} 

/ } 
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Appencfix  D:  Sample  SEU_CALC  output  file 


The  input  flies  used  in  this  run  were  those  listed  in  Appendix  B  and  C.  The  user 
selected  the  “low”  convergence  condition  by  entering  “1”  at  the  prompt.  The 
(HJtput  from  the  SEU_CALC  output  file  follows.  To  the  screen  the  program  writes 
some  adcfrtionai  intermecRate  results.  One  purpose  of  this  is  to  help  assure  the 
user  that  the  program  is  still  executing. 


SEO  rate  calculator  -  1-D  Integration  method 
Integral  over  experimental  cross-section  data. 
Version  1.0  D.L.  Chenette  8  November  1991 


Convergence  parameter  *  1.  Integration  consistency  “0.1 


LET  spectrum  used  in  this  run: 

LET  (MeV/mlcron)  Flux  (/m2  s  sr)  "LEASAT  1"  spectrum 


LET  (MeV/micron) 
l.OOOOe-02 
1.3388e-02 
1.7925e-02 
2.3998e-02 
3.2129e-02 
4.3015e-02 
5.7590e-02 
7.7103e-02 
1.0323e-01 
1.3820e-01 
1.8503e-01 
2.4773e-01 
3.3166e-01 
4.4404e-01 
5.9449e-01 
7.9592e-01 
1.0656e+00 
1.4267e+00 
1.9100e>00 
2.5572e+00 
3.4237e+00 
4.5837e+00 
6.1368e-i-00 
8.2161e4-00 
l.lOOOe+01 
LET  spectrum  set. 


Flux  (/m2  s  sr) 
2.1972e+01 
1.8184e+01 
1.4426e+01 
1.1121e+01 
8.0889e-t-00 
5.9511e+00 
4.3873e+00 
3.1657e+00 
2.2608e-»-00 
1.6468e+00 
1.1756e+00 
7.7379e-01 
4.0817e-01 
1.9998e-01 
1.0080e-01 
5.1504e-02 
2.6613e-02 
1.3688e-02 
6.8406e-03 
3.4166e-03 
1.5384e-03 
6.7240e-04 
1.3942e-04 
5.6336e-08 
O.OOOOe-^00 
25  values  tabulated 


Device  Upset  Cross-Section  Tabulation 

LET  (MeV/mg/cm2)  Cross-Section  (cm2/device)  */ 
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1^7  CfTOss-aection  (cm2/devlc«) 


l.OOQOe-01 

O.OOOOetOO 

S.OQOOe-Ol 

l.OOOOe-04 

i.OQflOe't^OO 

3.23006-03 

1 . leoQetoo 

4.34006-03 

a.OQOQe^ep 

5.7700e’-03 

i .  30OC!e-t^OG 

7. 20006^03 

S.GGQQe-^OQ 

8.7300e-03 

8.?GQ€te400 

l.OQOOe-02 

a.GQQOe-t^OO 

1.12006-02 

l.SOeOe^Ql 

1.1800e-02 

2 . lOQQe^Gl 

1.24Q0e-02 

3.00QOe401 

1.3000e-02 

3i.5GO0e-^Ol 

1.3900e-02 

4.2000e-^01 

1.4800e-02 

Sensitive  region  (microns) :  X  «  39.00 

Numl^er  of  bits  per 

device  -  10.24 

C^osa  section  set. 

14  values  tabulated 

1.00 


Integration  conqplete. 
Result  ">  ^.329Se-04 
Extrene  limit  bounds 
Result  *  9.5534e-0l 
Exbrei^  limit  bounds 


upsets  per  bit  per  day 
are  6.4874e-04  to  1.3410e-03 
upsets  per  device  per  day 
are  6.6431e-01  to  1.3732e+00 
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AppencNx  E:  Source  code  listings 


Listings  of  the  current  version  of  the  source  a>de  are  included  here.  All  of  the 
source  code  is  written  in  ANSI  C.  The  source  code  is  contained  in  three  files 
named  SEU_CALC.C,  PATHDtST.C.  and  INP_LIST.C.  All  of  the  source  code  is 
in  the  public  cksmain.  The  subroutines  contained  in  each  file  are: 

SEU_CALC.C 

maln() 

integrandO 

myfluxO 

IgaussO 

doJntO 

PATHDIST.C 

pIdsetO 

pWt() 

pld{) 

INP_LIST.C 

inputJistO 
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Page  1 


•EH  oac.c 

■oaSap,  Mareh  16,  1992  10:23  Ml 


finclude  <8tdio.h> 
fiaelude  <Mth.h> 

/*  Routines  to  do  ID  Integral  of  single-event  iqpset  cross-section 
over  pathlsffigth  distribution  and  LET  spectrum  to  calculate  the 
rate  of  8lngle-event-cg>sets. 

Written  1^:  D.L.  Cbenette  Version  1.3  8  Novonber  1991 

Required  source  code  modules:  SEUjCALC.C  INP_LIST.C  PATHDIST.C 

SEU  CALC:  this  module  with  main,  LET  spectnsn,  emd  cross  section  models 

nvjjST:  Infiut  module  to  read  DEVICE  and  LET_SPECTRUM  files 

PATHDIST:  module  to  calculate  pathlength  distribution  of  rectangular  solid 


To  compile  and  produce  an  executable  file  under  VMS,  do  the  following 
(make  sure  all  files  are  in  the  same  directory  before  you  start) : 
Define  the  pr^er  link  libraries  (may  be  done  in  the  login. com): 
DEriWE  LNK$LIBEVlRy  SYS$LIBRARy:VAXCCURSE.OLB 
EXEFIME  LNK8LZBRARy_l  SYS$LIBRARY:VAXCRTL.OLB 
Compile  all  modules  individually  with  the  following  commands  to 
generate  .OBJ  files  (unnecessary  if  .OBJ  file  exists  alreex^) : 

CC  PATBDIST 
CC  IWP  LIST 
CC  SEU  jaLC 

Mmo  this  Ts  done  you  should  have  .OBJ  files  of  each  of  these. 

Link  all  modules  to  make  the  executable  file  with  the  following: 

LIW  SEOjCALC,INP_LIST,  PATHDIST 
The  link  step  should  produce  the  executable  file  SEUjCALC.EXE 
which  you  can  run  with  the  command  RUH  SEUjCALC 

Before  RUHhlng,  however,  you  must  DEFINE  DEVICE  and  LET_SPECTRUM 
ii^nit  files  as  well  as  the  SEUjCALC  output  file.  Alternatively, 
these  file  name  may  be  used  directly  if  the  files  are  in  the 
current  directory. 

The  DEVICE  euKl  LET_SPECTRUM  input  files  are  nearly  free-fonnat. 

In  each  file  special  text  keys  are  used  to  distinguish  comments 
or  other  text  from  the  information  required  by  the  program. 

The  levies  file  describes  the  Integrated  circuit  part  for  the 
SlU  rate  calculation.  The  following  information  is  required: 

1.  Diswnsions  of  the  sensitive  cell  in  microns  (X,Y,Z). 

The  2  dlMnslon  Is  assumed  to  be  depth  into  the  silicon. 

The  quantity  x*Y  is  the  sensitive  surface  eurea  of  a  cell. 

2.  The  nunft>er  of  sensitive  volumes  (N)  in  the  cross-section 
table  which  follows.  If  the  cross-section  is  tabulated 
per  bit,  forexample,  N  should  be  set  to  l. 

3.  A  2-«>lumn  listing  of  I£T  and  upset  cross-section  valvms. 

The  mits  for  the  LET  must  be  NeV/micron,  and  the  units 
for  the  v^set  cross-section  must  be  in  square  centimeters 
per  device  (or  per  bit,  if  N  is  set  to  1,  see  above) 

The  keys  f<n:  X,  Y,  Z,  and  M  are  the  two-character  sets  X>  Y»  Z*  or 
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N“  The  first  numerical  value  found  after  the  key  on  the  same  line 
is  interpreted  as  the  value  for  the  required  parameter.  The  keys 
and  their  values  may  appe2u:  on  one  or  more  lines.  They  should 
appear  in  the  DEVICE  file  before  the  cross-section  key  text  string. 

For  the  cross-section  table  of  the  DEVICE  file  and  for  the  LET 
spectrijm  of  the  LET_SPECTRUM  file  the  key  is  the  3-characters  LET 
The  first  line  containing  "LET"  is  assumed  to  be  a  comment  and  is 
listed  to  the  output.  Subsequent  lines  to  the  end  of  the  file  are 
expected  to  contain  at  least  2  numbers  (and  may  contain  additional 
values  or  other  text) .  The  first  2  numbers  on  each  line  are  used 
and  interpreted  as  the  LET  and  the  upset  cross-section  (for  the 
DEVICE  file)  or  the  LET  and  the  integral  fl\ax  (per  square  meter 
per  second  per  steradian,  for  the  LET_SPECTRUM  file) .  If 
only  one  number  is  found  on  a  line,  that  line  will  be  igored  as 
a  comnent.  The  format  of  the  nuimbers  not  critical.  Fixed-point, 
floating  point  and  exponential  notations  are  recognized.  The  program 
lists  what  it  has  read  to  the  output  file.  When  using  a  new  input 
DEVICE  or  LET_SPECTRUM  file  this  listing  should  be  checked  to  make 
sure  that  it  is  being  interpreted  correctly. 

Restrictions:  The  first  entry  of  the  DEVICE  file  should  have  a 
non-zero  LET  with  a  cross-section  of  zero  to  start  the  inegration. 
The  both  files  must  be  monotonic  increasing  in  LET.  The  device  file 
must  be  monotonic  Increasing  in  cross-section  with  increasing  LET. 
The  LET_SPECTRUM  file  (Since  it  is  an  integral  flux)  must  be 
monotonic  decreasing  in  the  flux  with  increasing  LET. 

Meucimum  number  of  entries  in  upset  cross-section  is  now  set  to  50. 
Meutimum  number  of  entries  in  the  LET  spectrum  is  now  set  to  200. 

If  you  want  to  specify  a  step-function  upset  cross-section,  put  in 
a  two  line  cross-section  listing.  An  example  DEVICE  file  follows: 

X=  40  Y=40  Z=  1  Sensitive  volume  is  40x40x1  microns 

N”  1024blts  LET  for  step  function  example  for  IK  device 
1.0  0.0  step  function  upset  cross-section  with  critical 

1.0  0.001  LET  threshold  =  1  MeV/mg/sqcro  and  total  upset 

cross  section  for  the  device  of  0.001  sqcm 

User  input:  an  integer  (1,  2,  or  3)  used  to  control  the  convergence 
of  the  integration.  A  value  of  1  requires  convergence  to  10%,  2 
requires  convergence  to  1%,  and  3  requires  convergence  to  0.1%.  In 
most  cases  the  upset  cross-section  is  not  measured  so  precisely  that 
the  strictest  convergence  requirement  is  necessary  a  to  optimize  the 
accuracy  of  the  result.  Execution  is  fastest  using  1,  3  is  pretty 
slow  on  a  Mac  or  PC  but  is  probably  not  too  bad  on  a  mainframe. 

*/ 


/*  The  parameter  INTRP  is  the  number  of  Interpolation  steps  used  between 
listed  values  of  the  upset  cross-section.  If  faster  execution 
is  desired  this  may  lt)e  set  to  1.  If  a  finer  interpolation  of  the 
listed  cross-section  is  desired,  it  may  be  set  higher.  I  use  8. 

The  \q^r  and  lower  limits  of  the  upset  rate  axe  calculated  using 
the  obvious  upper  and  lower  bounds  of  the  upset  cross-section.  They 
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are  certainly  conservative  limits.  The  mid  value  Is  from  this 
Interpolation.  */ 
fdeflne  IlfERP  8 

double  Qc;  t*  critical  charge  In  PICOCOULOMBS 

tdeflne  NOMXST  50  /*  number  of  entries  In  following  table  *! 

double  xsli^ [2] [NOMXST ] ;  /*  cross-section  table  as  provided  */ 

fdeflne  NOMSPC  200  /*  number  of  LET  entries  In  table  below  */ 

doid>le  spectra [2] [NOMSPC] ;  /*  LET  spectrum  as  provided  */ 

double  Accuracy  0.001;  /*  End  the  Integration  when  two  successive 

Integrals  differ  by  less  than  this 
amount  as  a  fraction  of  the  value. 

Note  that  the  error  In  the  result  of  the 
Integral  likely  will  be  larger  than  this  */ 

I 

/*  The  following  constants  are  used  to  adjust  units  in  the  program  */ 
double  funlts  «•  0.0000000864;  /*  convert  from  /s  m2  to  /day  imi2  86400*le-12  */ 

double  Aotult  l.OeS;  /*  convert  area  In  squaure  cm  to  squaire  microns  */ 

double  imult  ~  0.233;  /*  convert  thickness  from  microns  to  mg/ cm2  (SI)  */ 

double  pCpMeV  «  0.0445;  {*  picoCoulombs  per  MeV  for  energy  deposit  in  Si  */ 

double  qnult;  /*  to  convert  E-  to  Q-dep.  MeV/mg/cm2  to  pC/micron. 

Defined  in  main  as  the  product  of  Tmult  *  pCpMeV  */ 

double  xsect [NOMXST] [2];  /*  xsinp  data  converted  to  my  units  */ 

double  spect [NOMSPC] [31 ;  /*  spectra  data  converted  to  my  units  */ 

double  X,  y,  z,  nblts,  Xs,  dQc,  dXs; 
double  rnt[INTRP],  ont[lNTRP]; 
int  nspect,  ndevxs; 

FILE  *po; 

malnO 

{ 

int  1,  j,  k,  ntrp,  index,  model; 

double  sqrt ( ) ,  pldset () ,  my f lux ( ) ; 

double  Imin,  Imaw,  eff,  lasteff,  lastleff,  lastueff; 

double  mid,  Imld,  upper,  lower,  result; 

char  kar,  title [NOMSPC] ; 

qpmlt  ■  linult  *  pCpMeV;  /*  define  conversion  from  MeV/mg/cm2  to  pC/um  */ 
po  -  fc^n("SEO_CALC","w"); 

printf("\f  SEO  rate  calculator  -  1-D  integration  method"); 
printf("\n  Integral  over  experimental  cross-section  data."); 
printf("\n  Version  1.0  D.L.  Chenette  8  November  1991  \n"); 

k  -  0; 
do  ( 

printf("\n\n  Enter  convergence  precision  parameter  (l«>low,  2»>ffledlum,  3*>high) : 
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scanf ("%d",&i) ; 
k  +=  1; 

)  while (k  <  4  fis  (i<l  | I  i>3) ) ; 
if  (k  >=  4)  { 

printf("\n  Input  error.  Got  %ci  Should  be  1  2  or  3"); 
goto  finish; 

) 

if  (i  »==  1)  Accuracy  =  0.1; 
if  (i  =  2)  Accuracy  =  0.01; 
if  (i  =»  3)  Accuracy  =  0.001; 

if  (po  !=  NULL)  { 

fprintf (po, "\f  SEU  rate  calculator  -  1-D  Integration  method"); 
f print f (po, "\n  Integral  over  experimental  cross-section  data."); 
fprintf (po, "\n  Version  1.0  D.L.  Chenette  8  November  1991  \n"); 
fprintf (po, "\n  Convergence  parameter  =  %d.  Integration  consistency  =  %g\n", 

1, Accuracy) ; 

} 

if  ((nspect  =  input_list ("LET_SPECTRUM", "LET", NUMSPC, spectra[0) , spectra[l] , title) )  >  1 
printf("\n  LET  spectrum  used  in  this  run:"); 
printf("\n  %s", title); 

printf("\n  LET  (MeV/mlcron)  Flux  (/m2  s  sr)"); 
if  (po  !“  NULL)  { 

fprintf (po, "\n  LET  spectrum  used  in  this  run:"); 
fprintf (po, "\n  %s",title); 

fprintf (po, "\n  LET  (MeV/micron)  Flux  (/m2  s  sr)"); 

) 

for (i«0;i<nspect;i++)  {  /*  convert  LET  spectrum  units  */ 

printf("\n  %12.4e  %12.4e", spectra[0) [i] , spectratl] [ij ) ; 

if  (po  !*  NULL)  { 

fprintf (po, "\n  %12.4e  %12.4e”, spectra[0) [i] , spectra[l] [i] ) ; 

} 

/"  convert  LET  (MeV/micron)  to  charge  dep  (pC/micron)  */ 
spect(i]  [0]  =  spectra [OJ  [i]  *  pCj^V; 
spect[l][l]  =  spectra[l] [i] ; 
if  (i  <  nspect-1)  spect[i][2]  = 

(spect[i+l] [1]  -  spectti] [l])/(specttl+l] [0]  -  spect [i] [0] ) ; 

/*  Note:  conversion  of  flux  to  per  sq  micron  and  day 
is  done  only  in  printout  at  end  to  keep 
numbers  larger  (comparable  to  1)  *! 
if  ((i  >  0)  &&  (spectra[0] [i-l]  >  spectra[0] [i]  II 
spectra [1 ] (i-lj  <  spectra [1] [i] ) )  { 
printf("\n  Error  in  LET_SPECTRUM  file.  Must  Ise  monotonic"); 
goto  finish; 

) 

} 

printf("\n  LET  spectrum  set.  %d  values  tabulated\n", nspect) ; 
if  (po  !»  NULL)  { 

fprintf (po, "\n  LET  spectrum  set.  %d  values  tabulatedXn", nspect) ; 

> 

) 
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alM  { 

printfCXn  LET  Spectnan  not  fovmd”); 
goto  finish; 

} 

If  ((lnput_llst ("DEVICE", "X-",0,tx, fix, title)  *=  -1)  fifi 
(input_li8t ("DEVICE", "Y-",0,fiy,fiy, title)  =  -1)  fifi 
(ir5>ut_llst ("DEVICE", "Z=",0,fiz,fiz, title)  =  -1)  fifi 
(ii^t_llst ("DEVICE",  "ll«",0,finbits,finbits, title)  -1)  &fi 

((ndevxs  -  input_list("DEVICE", "LET", NUMXST,xsinp[0],xsinp[l], title) )  >  0)  fifi 
X  >  0.0  fifi  y  >  0.0  fifi  z  >  0.0  fifi  nblts  >  0.0)  { 

printf("\n  Device  Upset  Cross-Section  Tribulation"); 
printf("\n  %s”, title); 

prlntf("\n  LET  (MeV/in9/cm2)  Cross-section  (aii2/device) ") ; 
if  (po  !-  NULL)  ( 

fprintf (po, "\n  Device  Upset  Cross-Section  Tabulation"); 
fprintf (po, "\n  %s", title); 

fprintf (po, "\n  LET  (MeV/mg/cmZ)  Cross-section  (cm2/devlce) ") ; 

} 

for(i«0;l<ndevxs;i++)  {  /*  convert  cross-section  data  */ 

prlntf("\n  %12.4e  %12.4e",xsinp[0] [i] ,xsinp[l] [i] ) ; 

if  (po  !-  NULL)  { 

fprintf (po,"\n  %12.4e  %12.4e",xsinp[0] [i] ,xsinp[lj [i] ) ; 

} 

/*  convert  LET  (MeV/mg/cm2)  to  charge  deposit  in  p>C  */ 
xsectli] [0]  »  xsinptO] (i]  "  <jBult  *  z; 

/*  cross-section/device  to  upset  efficiency  per  bit  */ 

/*  by  dividing  by  nianber  of  bits  and  area  per  bit  */ 
xsectCl] [1]  ■  xsinpCl] [i]  /  (x  "  y  *  nbits  /  Amult); 

/*  calculate  slopes  for  linear  interpolations  later  */ 
if  ((i  >  0)  fifi  (xsir^tO) [i-l]  >  xsinptO] [i]  I  I 
xsinptl]  [i-1]  >  xsir^tl]  tij ) )  ( 
printf("\n  Error  in  DEVICE  data  file.  Must  be  monotonic" ) ; 
goto  finish; 

} 

) 

printf("\n  Sensitive  region  (microns):  X  =  %.2f  Y  -  %.2f  Z  =  %.2f",x,y, z) ; 
printf("\n  Number  of  bits  per  device  =  %. Of”, nbits) ; 
printf("\n  Cross  section  set.  %d  values  tabulated\n\n", ndevxs) ; 
if  (po  !-  NULL)  { 

fprintf (po, "\n  Sensitive  region  (microns):  X=%.2f  Y-%.2f  Z**  %.2f",x,y, 

fprintf (po, "\n  Nvanber  of  bits  per  device  *  %. Of", nbits) ; 

fprintf (po, "\n  Cross  section  set.  %d  values  tabulated\n\n",nctevxs) ; 

) 

) 

else  ( 

prlntf("\n  Device  description  file  not  found  or  in  error"); 

prlntf("\n  Sensitive  region  (microns):  X  -  %.2f  Y  *  %.2f  Z  «  %.2f",x,y,2); 

prlntf("\n  Nunber  of  bits  per  device  «  %. Of”, nbits) ; 

prlntf(”\n  %d  cross-section  values  found  \n\n", ndevxs) ; 

goto  finish; 

} 
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t*  range  of  pathlength  Integral  is  0  to  maximum  pathlength  in  region  */ 
Iffiin  -  0.0; 

Imax  «  pldset(x,y, z) ; 

/* 

printf("\n  Press  return  to  continue"); 
scanf  (•%c",  s)car)  ; 

*/ 

lastleff  «  0.0; 
lastueff  =  0.0; 
lasteff  =  0.0; 
lower  »  0.0; 
upper  -  0.0; 
umld  >  0.0; 

Imld  -  0.0; 

for (i=0;i<ndevxs;i++)  {  /*  do  all  parts  of  cross-section  */ 

Qc  »  xsect[i] [01; 

Xs  =  xsect [i] [IJ ; 
if  (i  <  ndevxs-1)  { 

dQc  =  (xsect [i+1] [0]  -  xsect [i] [0] )/( (double)  INTRP) ; 
dXs  *  (xsect [i+1] [1]  -  xsect[i] [lJ)/( (double)  INTRP); 

} 

printf("  Part%3d  of%3d  Qc=%11.4e  Xs=%11.4e", 

i+1, ndevxs, Qc, Xs) ; 

)c  -  ((dQc  <=  0.0)  I  I  (i  =  ndevxs-1) )  ?  1  :  INTRP; 
for(  j*0;  j<)c;  j++)  {  /"do  interpolation  steps  */ 
model  *  do_int (Accuracy, Imin, Imax, rnt+  j , ont+  j ) ; 

Qc  +“  dQc; 

Xs  +«  dXs; 

) 

printf("  I»%11.3e\n",rnt[01 ) ; 

upper  +=  rnt[01  "  (xsect [i+1] [1]  -  lastueff); 

lastueff  =  xsect [i+1] [1] ; 

lower  +=  rnt[0]  *  (xsect[i][l}  -  lastleff); 

lastleff  =  xsect  [i][l]; 

Imld  +“  rnt[0]  *  (xsect[l] [1]  -  lasteff); 
lasteff  =  xsect [i]  [1]  +  dXs  *  ()c  -  1); 
if  ()t  >  1)  ( 

for  ( j=-0;  j<)c;  j++)  umld  +=  rnt[j]  *  dXs; 
for  ( j”l;  j<)c;  j++)  Imid  +=  rnt[j]  *  dXs; 

) 

) 

printf ("\n  L»%11.4e  U=%12.4e  lmid=%12.4e  uraicN=%12.4e", 
lower,  upper,  Imld,  umld) ; 
lower  *=  f units; 
upper  funits; 

umld  «  (umld  +  Imld)  *  funits  /  2; 
printf ("\n\n  Integration  cooqjlete.") ; 

printf ("\n  Result  -%12.4e  upsets  per  bit  per  day”,unld); 
printf ("\n  Extrone  limit  iDounds  are%12.4e  to%12. 4e", lower, upper) ; 
printf ("\n  Result  »%12.4e  upsets  per  device  per  day", vBnid*nbits) ; 
printf ("Xn  Extreme  limit  ]x)unds  are%12.4e  to%12.4e", 
lower "nbits, upper"nblts) ; 
if  (po  !-  NULL)  ( 
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fprintf (po, "\n\n  Integration  conqplete.”) ; 
fprintf  (po, ”\n  Restilt  ■•%12.4e  upsets  per  bit  per  day",uRild); 
fprintf (po, ”\n  Extreme  limit  bounds  are%12.4e  to%12. 4e”, lower, upper) ; 
fprintf (po, ”\n  Result  »%12.4e  upsets  per  device  per  day",utnld*nblts); 
fprintf (po, ”\n  Extreme  limit  boxmds  are%12.4e  to%12.4e”, 
lower*nblts,u{^r*nblts) ; 

} 

finish: 

printf("\n  End  of  run.  \n\n"); 
if  (po  !»  NULL)  fclose(po); 

} 

double  integrand(l)  /*  this  is  the  function  called  by  twodlntO  */ 
double  1; 

( 

double  pldt ( ) ,  myf lux ( ) ; 

If  (1  <-  0.0)  ret\im(0.0); 
return (pldt (1)  *  myflux(Qc/l) ) ; 

> 

double  myflux(L) 
double  L; 

{ 

Int  i; 

static  int  old!  »  0; 

if  (L  >-  specttoldl] [0]  L  <*  spect[oldi+ll [0] )  { 

1  -  oldi; 

> 

else  { 

if  (L  <-  0.0  I  I  L  >”  spect [nspect-1] [0] )  return(O.O); 
for (l*n8pect-2;l>0;i — )  if  (L  >  spect(ij[0))  break; 

) 

oldi  *1; 

return (spect [i] [1]  +  spect[l][2J  *  <L  -  spect [1 J [OJ )) ;  /*  calc  flux  at  L  */ 

) 

dovible  lgauss(a,b)  /*  do  Gaussian  Integral  of  integrand  from  a  to  b  */ 

double  a,b; 

( 

int  j; 

double  xr,xm,dx,s; 
double  IntegrandO; 

static  double  x[]-(0. 0,0. 1488743389,0. 4333953941, 

0 . 6794095682, 0 . 8650633666, 0 . 97390652} ; 
static  double  wt]-( 0.0, 0.2955242247, 0.2692667193, 

0 . 2190863625, 0 . 1494513491, 0 . 06667134 } ; 

xa)>0.5*(bfa)  ; 
xr«0.5f  (b-a)  ; 


/*  Calculate  flux  at  LET  =  L  */ 
/*  L  in  units  of  pC/micron  */ 
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for  (j-1; j<-5; j++)  { 
dx-9tr*x[  j]  ; 

s  +■  wt  jj  *  (integrand(wn+dx)  +integrand(xm-dx) ) ; 

) 

return (s  xr) ; 

} 

Int  do_lnt  (AcCfRtn, mx,  res,  old) 
double  Acc,  inn,  rax,  *res,  *old; 

{ 

dovible  dx,  ans%fez,  panswer,  xmln,  xraax; 
double  min,  max,  eps,  ddx,  trlnt,  ntrlnt; 
double  rpart,  opart,  IntegrandO,  IgaussO; 
int  i,  section,  ntr; 

ntr  »  20; 

ntrint  ■  (double)  ntr; 
h 

printf("\n  Integrate  from  %12.4e  to  %12.4e\n",mn,mx) ; 
if  (po  !-  NULL)  { 

fprintf (po, "Xn  Integrate  from  %12.4e  to  %12.4e\n”,mn,mx) ; 

) 

*( 

dx  «  (rax  “  ran)  /  ntrint; 
for(i«0;l<-ntr;l++)  ( 
xmln  «  xraax; 

xraax  -  ran  dx  *  ((double)  i); 
ddx  •  eps; 

eps  *  1 ntegr and (xraax ) ; 
if  (1  —  0)  { 
trint  =  0.0; 

) 

else  ( 

trint  +=  dx  *  (ddx  +  eps)  /  2.0; 

) 

/* 

printf("  X  =%12.4e  f  »%12. 4e\n", xraax, eps) ; 
if  (po  !“  NULL)  { 

fprintf(po, "  X  «%12.4e  f  «%12. 4e\n", xraax, eps) ; 

) 

*/ 

) 

/* 

printf("  %d  point  trapedoid  Integral  »%12. 4e\n", ntr, trint ) ; 
if  (po  !-  NULL)  ( 

fprintf(po, "  %d  point  trapedoid  integral  “%12. 4e\n", ntr, trint) ; 
) 

*/ 

rpart  *  0.0; 

opart  -  0.0; 

answer  »  0.0; 

ddx  -  (nx  -  ran)  /  ntrint; 

for  ( section^;  8ectlon<ntr ;  secti<m-t"f )  { 
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/* 

*/ 


*/ 

I* 

*/ 


■in  »  an  +  *  ( (doiA>le)  section) ; 

met  -  aii^  +  d^* 
foraL-4,;i<<$00;i^)  { 

•  (pen  ~  Biin)  /  ((double)  1); 

Rpnaiier  «  ansMer; 
anawiMc  -  0.0; 
xp^  aIa; 

%dille(xain  <  max)  { 
xpax  «  xpln  +  dx; 

1#  (xpaa  >  max)  xmax  •«  max; 
pnamar  ^  igau8s(xmln,xpax); 
xpin  «  xmax; 

) 

If  ((1  >  1  fit  answer  !«  0.0)  II  (1  >  2))  ( 

if  (answer  !»  0.0)  eps  -  (answer  -  panswer)  /  emsw.er; 
else  affs  «  -1.0; 

^intfC  %5d  %12.4e  %12.4e  %])2. 4e\n‘'>if]^swer, answer, eps) ; 
if  (po  !«  IKJLL)  { 

fprlntf(po,”  %5d  %12.4e  %12.4e  %12.4e\n*, 

1,  panswer,  ans%ier,  eps) ; 

) 

if  (fabs(eps)  <  Acc  II  fabs(eps*answer/trlnt)  <  Acc/10.0)  break; 
) 

) 

xpart  *•  answer; 
opart  >«•  panswer; 

) 

*old  *•  qpart; 

*rea  «  xpart; 
return (i/2) ; 
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/*  INPOT_LIST.C 

Purpose:  Read  any  listings 

Coded  by  O.L.  Chenette,  8  November  1991 

Change  log: 

*/ 


ilnclude  <stdlo.h> 
ilnclude  <stdllb.h> 
tinclude  <Bath.h> 

idefine  MAXLLEN  250 

int  lnput_list (if name, key,maxn, coll, col2, title) 
char  *lfnaine,  *key,  *title; 
int  maxn; 

double  *coll,  *col2; 

{ 

char  s [MAXLLEN] ,  *c,  *d; 
int  k; 

double  *cl,  *c2; 

FILE  *fp; 

k  -  0; 

if  (ifname  (fp  »  fopen(ifname, "r") )  !*  MOLL)  { 

/* 

prlntf("\n  Reading  data  from  file:  %s”, ifname); 

*/ 

if  (*key)  { 

while (fgets(s, MAXLLEN, fp))  ( 
c  -  s; 
d  «  key; 

while (*c  &&  *c  !■  *d)  C++; 
while (*d  *c  it  (*c  “  *d))  { 
c  +=  1; 
d  +-  1; 

) 

if  (maxn  —  0  4t  *d  —  '\0’)  { 

while(*c  44  (*c  <  '-•  II  *c  >  '9')) 
if  (*c)  *coll  -  atof(c); 
return (-1) ; 

) 

if  (*d  “  '\0')  ( 
c  «  title; 
d  -  s; 
k  -  80; 

while ((*c++  -  *d++)  44  --k); 

*c  “  '\0'; 
k  -  0; 
goto  begin; 

) 


Version  1.0 


C++; 
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} 

priUtf  ("Vn  Input  fll«  %•  did  not  contain  kay  string  %s”,ifnaM,tey) ; 
goto  done; 

) 

begin: 

cl  >■  ooll; 
c2  ■■  ooI2; 

while  ffgets(s,MaLLEii,fp))  ( 

c  -  s; 

while (*c  «c  (*c  <  •-•  II  >  ‘g*))  C++; 
if  (*c)  { 

*cl  ■  atof(c); 

while  <*c  cc  (*c  >-  •-•  II  *c  —  ’+•))  C++; 
tdille(*c  sc  (*c  <  II  *c  >  '9'))  C++; 
if  (*c>  ( 

*e2  -  atof(c); 

Cl  +-  1; 
c2  +-  1; 

if  (++k  <—  naxn)  goto  done; 

i 

) 

) 

) 

done: 

if  (fp  *—  NOLL) 

prlRtf(”\n  Iig>ut  file  naned  %s  cannot  be  opened  for  reading”, ifname) ; 
else  fclose (fp) ; 

return (k); 

} 
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flnclude  <stcllo.h> 
flnclude  <8tdlib.h> 
tlnclvide  <math.h> 

dotible  cllst(1002];  /*  pathlength  distribution  calculator  method  taken  from  */ 
double  snax;  /*  Pickel  &  Blandford,  IEEE  Trans. Nuc.Sci.NS-27,  Apr. 1980  */ 

/*  for  right  rectangular  prism  geometry  */ 

/*  coded  by  D.  L.  Chenette  25  June  1987  */ 

double  pldset(zzl,  zz2,  zz3)  /*  set  up  pathlength  distribution  table  *! 
double  zzl,  zz2,  zz3;  /*  PLO  is  tabulated  every  0.1%  of  maximum  */ 

{  /*  pathlength  (1001  points  including  meuc)  */ 

int  i;  /*  returns  maxlmisn  pathlength  in  volvune  */ 

double  pldO,  sqrtO; 

smax  “  sqrt(zzl  *  zzl  +  zz2  *  zz2  +  zz3  *  zz3); 
for(l-0;i<1001;i++)  { 

dist[l]  -  pld( (double)  smax*i/1000.0,  zzl,  zz2,  zz3) ; 

) 

distdOOll  -  0.0; 

prlntf("\n  Pathlength  distribution  set  for  right  rectangular  prism."); 
printf("\n  Dimensions  of  sensitive  region:  x  =  %g  y  =  %g  z  =  %g", 
zzl,  zz2, zz3) ; 

printf("\n  Maximvim  pathlength  through  region  *  %g  \n\n”,smax>; 
return (smax) ; 

) 

double  pldt(s)  f*  return  tabulated  value  of  pathlength  distribution  */ 
double  s;  I*  for  pathlength  s  through  right  rectangular  prism  */ 

{  /*  must  have  called  pldset(x,y, z)  first  to  setup  table  */ 

int  i; 
double  x; 

if  (s  <  0.0  II  s  >  smax)  return(O.O); 

X  -  1000.0  *  s  /  smax; 
i  -  (int)  x; 

return ( dl St [i]  +  (dist[i+l]  -  dlst[l])  *  (x  -  i)); 

) 

double  pld(s, zzl, zz2, zz3)  /*  return  pathlength  distribution  for  rect  solid  *! 

double  s, zzl, zz2, zz3;  /*  s  »  pathlength,  zzl,  zz2,  zz3  are  dimensions  *i 

{  I*  pld  is  normalized  s.t.  Integral  over  all  s  is  *! 

int  i;  /*  equal  to  G-factor  -  4*pl* (average  proj.  area)  */ 

double  atanO ,acos() ,sqrt 0 ,pow() ;  /*  G-f  »  2*pi* (zzl*zz2+zzl*zz3+zz2*zz3)  */ 

double  X,  y,  z,  r,  r2,  g,  qsz,  fs; 

double  x2,  y2,  z2,  xy2,  xz2,  xy,  xz,  s2,  s3,  ssz,  ssz2; 
fs  •  0.0; 

if  (zzl  <■  0.0  II  zz2  <*0.0  II  zz3  <«  0.0  II  s  <  0.0)  return (fs); 
if  (s  >  sqrt(zzl*zzl  +  zz2*zz2  *  zz3*zz3) )  return(fs); 

s2  «  8  *  s; 

83  -  8  *  82; 

for (i«0;l<6;i++)  { 
if  (i  —  0)  { 
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X  »  zzl; 
y  -  zz2; 
z  «  zz3; 

) 

•Ise  if  (i  —  1)  ( 

X  «  zzl; 
y  •  zz3; 
z  »  zz2; 

else  If  (1—2)  { 

X  -  zz2; 
y  -  zzl; 
z  -  zz3; 

) 

else  if  (i  —  3)  { 

X  -  zz2; 
y  -  zz3; 
z  *  zzl; 

) 

else  if  (i  —  4)  { 

X  »  zz3; 
y  -  zzl; 
z  -  zz2; 

) 

else  if  (1  -=  5)  { 

X  -  zz3; 
y  »  zz2; 
z  «  zzl; 

> 

x2  »»  X  *  x; 

y2  -  y  *  y; 

z2  *  z  *  z; 

r2  =  x2  +  y2  +  z2; 

r  «  sqrt (r2) ; 
xz2  “  x2  +  22; 
xz  *  sqrt(xz2); 
xy2  *  x2  +  y2; 
xy  *  sqrt(xy2); 
if  (s  >  z)  { 

ssz2  *  s2  -  z2; 
ssz  -  3qrt(ssz2); 

) 

else  { 

ssz2  »  0.0; 
ssz  »  0.0; 

) 

if  (s2  >  xz2)  qsz  *  sqrt(s2  -  xz2} : 
else  qsz  -  0.0; 
if  (s  <  z)  { 

g  -  y2/(xz2  *  r)  +  2.0  *  (r  -  xz)/(xz  *  r)  +  r/xy2  -  xz/x2 
+  y2  *  z/(x2  *  xy2)  -  1.125  *  y2  *  s/ (r2  *  xz2); 
g  *-  2.0  *  x2  /  3.0; 


) 
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else  if  (s  <  xz)  ( 

g  »  y2/(xz2  *  r)  +  2.0  ♦  <r  -  xz)/{xz  *  r)  +  r/xy2  -  xz/x2 

-  1.125  *  y2  *  s/(r2  *  xz2) 

-  ssz  *  ssz2  *  (x  +  (y2  -  x2)/xy)  /  (s3  *  x2) 

+3.  *  z2  *  y2  *  (1.  -  z2/{3.  *  s2))/{2.  *  s  *  x2  *  xy2) 

+  1.125  *  ssz2  *  ssz2  *  y2  /  (x2  *  s3  *  xy2) ; 
g  *-  2.0  *  x2  /  3.0; 
g  +-  X  *  y  *  z2  *  atem(y/x)  /  s3 

-  z2  *  ssz  *  ((y2  -  x2)/xy  +  x)/s3  +  z2*ssz2*y2/ (2.*s3*xy2) ; 

) 

else  if  (s  <  r)  ( 

g  -  y2/(xz2  *  z)  -  y  *  qsz/(s  *  xz2)  -  2./r  +  r/xy2 

-  s/ssz2  +  2./s  -  1.125*s* (xy2/r2  -  ssz2/s2)/z2 

-  ssz  *  ssz2  *  ((y2-x2)/xy  -  (y*qsz  -  x2) /ssz) / (s3*x2) 
+(3./(2.*x2))*(0.75*s  -  z2/(2.*s)  +  5.*z2*z2/ (12.*s3) ) 

* (y2/xy2  -  (s2  -  xz2)/ssz2); 

g  2.0  *  x2  /  3.0; 

g  +-  (x*y*z2/ (s*s2) ) * (atan(y/x)  -  acos(x/ssz)) 

-  (z2/a3) *SS2* ( (y2-x2) /xy  -  (y*qsz  -  x2)/ssz) 

+  (z2/(2.*s3))  *  ssz2  *  (y2/xy2  -  (s2-xz2) /ssz2) ; 

) 

else  g  >  0.0; 
fs  +-  g; 

) 

return (8.0  *  fs) ; 

) 
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